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Abstract 

Noncommutative geometry has been slowly emerging as a new paradigm of 
geometry which starts from quantum mechanics. One of its key features is 
that the new geometry is spectral in agreement with the physical way of 
measuring distances. In this paper we present a detailed introduction with 
an overview on the study of the quantum nature of space-time using the 
tools of noncommutative geometry. In particular we examine the suitability 
of using the spectral action as action functional for the theory. To demon- 
strate how the spectral action encodes the dynamics of gravity we examine 
the accuracy of the approximation of the spectral action by its asymptotic 
expansion in the case of the round sphere S 3 . We find that the two terms 
corresponding to the cosmological constant and the scalar curvature term 
already give the full result with remarkable accuracy. This is then applied 
to the physically relevant case of S' 3 x S 1 where we show that the spectral 
action in this case is also given, for any test function, by the sum of two 
terms up to an astronomically small correction, and in particular all higher 
order terms a<m vanish. This result is confirmed by evaluating the spectral 
action using the heat kernel expansion where we check that the higher order 
terms 04 and a% both vanish due to remarkable cancelations. We also show 
that the Higgs potential appears as an exact perturbation when the test 
function used is a smooth cutoff function. 
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1. An overview 

Our experimental information on the nature of space-time is based on two 
sources: 

• High energy physics based on cosmic ray information and particle 
accelerator experiments, whose results are encapsulated in the Stan- 
dard Model of particle physics. 

• Cosmology based on astronomical observations. 

The large scale global picture is well described in terms of Riemannian ge- 
ometry and general relativity, but this picture breaks down at high energy 
where the quantum effects take over. It is thus natural to look for a par- 
adigm of geometry which starts from the quantum framework, where the 
role of real variables is played by self-adjoint operators in Hilbert space. 
Such a framework for geometry has been slowly emerging under the name 
of noncommutative geometry. One of its key features, besides the ability 
to handle spaces for which coordinates no longer commute with each other, 
is that this new geometry is spectral. This is in agreement with physics in 
which most of the data we have, either about the far distant parts of the 
universe or about high energy physics, are also of spectral nature. The red 
shifted spectra of distant galaxies or the momentum eigenstates of outgo- 
ing particles in high energy experiments both point towards a prevalence of 
spectral information. In the same vein the existing unit of time (length) is 
also of spectral nature. From the mathematical standpoint it takes some 
doing to obtain a purely spectral (Hilbert space theoretical) counterpart of 
Riemannian geometry. One reason for the difficulty of this task is that, as 
is well known since the examples of J. Milnor [I], non- isometric Riemannian 
spaces exist which have the same spectra (for the Dirac or Laplacian oper- 
ators). Another reason is that the conditions for a (compact) space to be a 
smooth manifold are given in terms of the local charts, whose existence and 
compatibility is assumed, but whose intrinsic meaning is more elusive. 
The paradigm of noncommutative geometry is that of spectral triple. As 
its name indicates it is of spectral nature. By definition a spectral triple is 
a unitary Hilbert space representation of "something". This something is 
an equipment that allows one to manipulate algebraically coordinates and 
to measure distances. The algebra of the coordinates is denoted by A and 
is an involutive algebra, with involution a i— ► a*. The equipment needed 
to measure distances is the inverse line element D which is unbounded and 
fulfills D = D* . Altogether these data fulfill some algebraic relations, e.g. 
if we talk about the simplest geometric space i.e. the circle S 1 the relation 
between the complex unitary coordinate U and the inverse line element D 
is just [D, U] = U, which is in the vein of the Heisenberg commutation 
relations. 

Thus, a geometry is given as a Hilbert space representation of the pair (.4, D) 
and can be encoded by the spectral triple (A, H, D) where TL is the Hilbert 
space in which both the algebra A and the inverse line element D are now 
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concretely represented, the latter as an unbounded self-adjoint operator. 
This picture shares with the Wigner paradigm for a particle as an (irre- 
ducible) representation of the Poincare group the feature that it separates 
the kinematical relations from the choice of the Hilbert space representation. 
It is only when the latter is chosen that actual measurements of distances 
between points x and y can be performed by formulas such as 

Distance (x,y) = sup \f(x)-f(y)\, f G A, ||[A/]||<1 

where indeed the norm /]|| is the operator norm in Hilbert space and 
depends on the specific choice of the representation. 

We now have at our disposal a reconstruction theorem (c/ [2]) which shows 
that ordinary Riemannian spaces are neatly characterized among spectral 
triples by the following kinematical relations: 

• The algebra A is commutative. 

• The commutator [[D, a],b] = for any a J) £ A. 

• The following "Heisenberg type" relatio rfl holds] , for some a" £ A: 

(i) ^^[[Ao?],[A^],---,[A<]] = i 

a 

together with the following spectral requirements: 

• The k-th characteristic value of the resolvent of D is 0{k~ l l n ). 

• Regularity. 

• Absolute continuity. 

We refer to [2] for the precise statement. The meaning of ([I]) is that the 
determinant of the metric g^ v does not vanish, and more precisely that its 
square root multiplied by the volume form Y^ a a^daf A daf A • • • da® gives 1. 
The reason for the last two spectral requirements is technical and allows one 
to specify the regularity (C°°, real analytic.) of the space and to control 
the spectral measures. The first of the spectral requirements is crucial in 
that it bounds the "effective dimension" of the spectrum of the space in the 
representation. There are good physics reasons to consider that the appar- 
ent dimension, equal to four, of space-time is governed by the asymptotic 
behavior of the eigenvalues of the line element, which is the Euclidean prop- 
agator. Moreover this spectral dimension is not restricted to be an integer 
a priori and can model fractal dimension easily. The above reconstruction 
Theorem shows furthermore that the operator D in the spectral triple is 
a Dirac type operator, i.e. an order one operator with symbol given by a 
representation of the Clifford algebra. The restriction to spin manifolds is 
obtained by requiring a real structure i.e. an antilinear unitary operator J 

^Here the multiple commutator is defined as 

[Ti,T 2 , . . . ,T n ] = 2^ T CT (i)T CT (2) ■ ■ ■ T CT(n ) 

2 We assume for simplicity that the dimension n is odd 
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acting in Ti which plays the same role and has the same algebraic proper- 
ties as the charge conjugation operator in physics. When the dimension n 
involved in the reconstruction Theorem is even (rather than odd) the right 
hand side of ([IJ) is now replaced by the chirality operator 7 which is just a 
Z/2-grading in mathematical terms. It fulfills the rules 

(2) 7 2 = 1 , [7, a] = 0, a G A 
The following further relations hold for D, J and 7 

(3) J 2 = s , DJ = e'JD, J 7 = £"7 J, D7 = -jD 

where £,e',e" £ {—1, 1}. The values of the three signs e,e' ,e" depend only, 
in the classical case of spin manifolds, upon the value of the dimension n 
modulo 8 and are given in the following table [3]: 
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-1 


-1 


-1 
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-1 
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-1 
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e" 


1 




-1 




1 




-1 





In the classical case of spin manifolds there is thus a relation between the 
metric (or spectral) dimension given by the rate of growth of the spectrum 
of D and the integer modulo 8 which appears in the above table. For more 
general spaces however the two notions of dimension (the dimension modulo 
8 is called the iTO-dimension because of its origin in iT-theory) become 
independent since there are spaces F of metric dimension but of arbitrary 
ifO-dimension. More precisely, starting with an ordinary spin geometry M 
of dimension n and taking the product M x F, one obtains a space whose 
metric dimension is still n but whose i^O-dimension is the sum of n with 
the iTO-dimension of F, which as explained can take any value modulo 8. 
Thus, one now has the freedom to shift the iTO-dimension at very little 
expense i.e. in a way which does not alter the plain metric dimension. As 
it turns out the Standard Model with neutrino mixing favors the shift of 
dimension from the 4 of our familiar space-time picture to 10 = 4 + 6 = 2 
modulo 8 [I], [5]. The shift from 4 to 10 is a recurrent idea in string theory 
compactifications, where the 6 is the dimension of the Calabi-Yau manifold 
used to "compactify" . Effectively the dimension 10 is related to the existence 
of Majorana-Weyl fermions. The difference between this approach and ours 
is that, in the string compactifications, the metric dimension of the full 
space-time is now 10 which can only be reconciled with what we experience 
by requiring that the Calabi-Yau fiber remains unnaturally small. In order 
to learn how to perform the above shift of dimension using a 0-dimensional 
space F, it is important to classify such spaces. This was done in [6], [7j. 
There, we classified the finite spaces F of given ETO-dimension. A space F 
is finite when the algebra Af of coordinates on F is finite dimensional. We 
no longer require that this algebra is commutative. The first key advantage 
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of dropping the commutativity can be seen in the simplest case where the 
finite space F is given by 

(4) A = M k (C), H = M k (C), D = 0, JZ = C, £ G TLf 

where the algebra A = M k (C) is acting by left multiplication in TL = M k (C). 
We have shown in [8] that the study of pure gravity on the space M x F 
yields Einstein gravity on M minimally coupled with Yang-Mills theory for 
the gauge group SU(&). The Yang-Mills gauge potential appears as the inner 
part of the metric, in the same way as the group of gauge transformations 
(for the gauge group SU(fc)) appears as the group of inner diffeomorphisms. 
One can see in this Einstein- Yang-Mills example that the finite geometry 
fulfills a nice substitute of commutativity (of A) namely 

(5) [a,6°] = 0, Va,beA 

where for any operator a in 7i, a = Ja*J _1 . This is called the order 
zero condition. Moreover the representation of A and J in TL is irreducible. 
This example is (taking 7 = 1) of i<CO-dimension equal to 0. In [6] we 
classified the irreducible (A, TL, J) and found out that the solutions fall into 
two classes. Let Ac be the complex linear space generated by A in C(Tt), 
the algebra of operators in TL. By construction Ac is a complex algebra and 
one only has two cases: 

(1) The center Z (Ac) is C, in which case Ac = M k (C) for some k. 

(2) The center Z (Ac) is C C and Ac = M k (C) © M k (C) for some k. 

Moreover the knowledge of Ac = M k (C) shows that A is either M k (C) 
(unitary case), M k (M.) (real case) or, when k = 21 is even, M^(H), where 
EI is the field of quaternions (symplectic case). This first case is a minor 
variant of the Einstein- Yang-Mills case described above. It turns out by 
studying their Z/2 gradings 7, that these cases are incompatible with KO- 
dimension 6 which is only possible in case (2). If one assumes that one is in 
the "symplectic-unitary" case and that the grading is given by a grading of 
the vector space over H, one can show that the dimension of TL which is 2k 2 
in case (2) is at least 2 x 16 while the simplest solution is given by the algebra 
A = M2(H) © M^(C). This is an important variant of the Einstein- Yang- 
Mills case because, as the center Z (Ac) is C © C, the product of this finite 
geometry F by a manifold M appears, from the commutative standpoint, 
as two distinct copies of M. We showed in (()| that requiring that these two 
copies of M stay a finite distance apart reduces the symmetries from the 
group SU(2) x SU(2) x SU(4) of inner automorphismsj to the symmetries 
U(l) x SU(2) x SU(3) of the Standard Model. This reduction of the gauge 
symmetry occurs because of the second kinematical condition [[D,a],&] = 
which in the general case becomes: 

(6) [[D,a],b°] =0, Va,b£A 



'of the even part of the algebra 
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Figure 1 . Cutoff function f 

Thus the noncommutative space singles out 4 2 = 16 as the number of phys- 
ical fermions, the symmetries of the standard model emerge, and moreover, 
as shown in [9], the model predicts the existence of right-handed neutrinos, 
as well as the see-saw mechanism. In the above Einstein- Yang-Mills case, 
the Yang-Mills fields appeared as the inner part of the metric in the same 
way as the group of gauge transformations (for the gauge group SU(/c)) ap- 
peared as the group of inner diffeomorphisms. But in that case all fields 
remained massless. It is the existence of a non-zero D for the finite space F 
that generates the Higgs fields and the masses of the Fermions and the W 
and Z fields through the Higgs mechanism. The new fields are computed 
from the kinematics but the action functional, the spectral action, uses in 
a crucial manner the representation in Hilbert space. In order to explain 
the conceptual meaning of this spectral action functional it is important 
to understand in which way it encodes gravity in the commutative case. 
As explained above the spectrum of the Dirac operator (or similarly of the 
Laplacian) does not suffice to encode an ordinary Riemannian geometry. 
However the Einstein-Hilbert action functional, given by the integral of the 
scalar curvature multiplied by the volume form, appears from the heat ex- 
pansion of the Dirac operator. More generally it appears as the coefficient 
of A 2 in the asymptotic expansion for large A of the trace 

(7) Tt(f(D/A)) ~ 2A 4 / 4 ao + 2A 2 / 2 a 2 + / a 4 + . . . + A~ 2k f-2k^ + 2k + ■■■ 

when the Riemannian geometry M is of dimension 4, and where / is 
a smooth even function with fast decay at infinity. The choice of the 
function / only enters in the multiplicative factors fa = / °° f(u)u 3 du, 
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h = / °° f(u)udu, f = /(O) and f_ 2k = (-l) k « /(*)(0), i.e. the deriva- 



tives of even order at 0, for k > 0. Thus, when / is a "cutoff' function 
(c/ Figure [1]) it has vanishing Taylor expansion at and the asymptotic 
expansion ([7]) only has three terms: 



The term in A 4 is a cosmological term, the term in A 2 is the Einstein-Hilbert 
action functional, and the constant term CI4 gives the integral over M of cur- 
vature invariants such as the square of the Weyl curvature and topological 
terms such as the Gauss-Bonnet, with numerical coefficients of order one. It 
is thus natural to take the expression Tr(/(D/A)) as a natural spectral for- 
mulation of gravity. We are working in the Euclidean formulation i.e. with a 
signature (+, +, +, +) and the Euclidean space-time manifold is taken to be 
compact for simplicity. In the non-compact case we have shown in [10] how 
to replace the simple counting of eigenvalues of \D\ of size < A giver@ by 
(|7j), by a localized counting. This simply introduces a dilaton field. We also 
tested this idea of taking the expression Tr(/(D/A)) as a natural spectral 
formulation of gravity by computing this expression in the case of manifolds 
with boundary and we found [11] that it reproduces exactly the Hawking- 
Gibbons [12J additional boundary terms which they introduced in order to 
restore consistency and obtain Einstein equations as the equations of motion 
in the case of manifolds with boundary. Further, Ashtekar et al [13] have 
recently shown that the use of the Dirac operator in a first order formalism, 
which is natural in the noncommutative setting, avoids the tuning and sub- 
traction of a constant term. One may be worried by the large cosmological 
term A /4CI4 that appears in the spectral action. It is large because the value 
of the cutoff scale A is dictated, roughly speaking, by the Planck scale since 
the term A 2 /2fl2 is the gravitational action 16 ^. G / R^fgd^x. Thus it seems 
at first sight that the huge cosmological term A 4 /4CI4 overrides the more sub- 
tle Einstein term A 2 /2a2- There is, however, and even at the classical level 
to which the present discussion applies a simple manner to overcome this 
difficulty. Indeed the kinematical relation ([!]) in fact fixes the Riemannian 



Thus, if we vary the metric with this constraint we are in the context of 
unimodular gravity [J3], and the cosmological term cancels out in the com- 
putation of the conditional probability of a gravitational configuration with 
total volume V held fixed. The remaining unknown, then, is the distribution 
of volumes dfi(V), which is just a distribution on the half-line M+ 3 V . The 
striking conceptual advantages of the spectral action are 



for / a cutoff function 
'up to a numerical factor 



(8) 



Tr(/(ZVA)) ~2A 4 / 4 a + 2A 2 



/ 2 a 2 + /(0)a4 




a 
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• Simplicity: when / is a cutoff function, the spectral action is just 
counting the number of eigenstates of D in the range [—A, A]. 

• Positivity: when / > (which is the case for a cutoff function) the 
action Tr(/(Z)/A)) > has the correct sign for a Euclidean action. 

• Invariance: one is used to the diffeomorphism invariance of the grav- 
itational action but the functional Tr(/(D/A)) has a much stronger 
invariance group, the unitary group of the Hilbert space TC. 

One price to pay is that, as such, the action functional TY(/(-D/A)) is not 
local. It only becomes so when it is replaced by the asymptotic expansion 
([8j). This suggests that one should at least compute the next term in the 
asymptotic expansion (even though this term appears multiplied by the 
second derivative /"(0) = when / is a cutoff function) just to get some 
idea of the size of the remainder. In fact both D and A have the physical 
dimension of a mass, and there is no absolute scale on which they can be 
measured. The ratio D/A is dimensionless and the dimensionless number 
that governs the quality of the approximation (jSj) can be chosen to just 
be the number N(A) of eigenvalues A of D whose size is less than A, i.e. 
|A| < A. When / is a cutoff function the size of the error term in (jSj) should 
be 0(N~ k ) for any positive k, using the flatness of the Taylor expansion of 
/ at 0. In the case of interest, where M is the Euclidean space-time, a rough 
estimate of the size of N is the 4-dimensional volume of M in Planck units 
i.e. an order of magnitud^l of N ~ 10 214 (at the present radius, see section 
two for details). Thus, even without the vanishing of /"(0), the rough error 
term N^ 1 / 2 ~ 10~ 107 is quite small in the approximation of the spectral 
action by its local version ([Sj) . We shall in fact show that a much better 
estimate holds in the simplified model of Euclidean space-time given by the 
product S^x Sp. Another advantage of the above spectral description of the 
gravitational action is that one can now use the same action Tr(/(D/A)) 
for spaces which are not Riemannian. The simplest case is the product 
of a Riemannian geometry M (of dimension 4) by the finite space F of 
([4j). The only new term that appears is the Yang-Mills action functional of 
the SU(A;) gauge fields which form the inner part of the metric. This new 
term appears as an additional term in the coefficient 04 of A , and with 
the positive sign. In other words gravity on the slightly noncommutative 
space M x F gives ordinary gravity minimally coupled with SU(fe)-Yang- 
Mills gauge theory. The latter theory is massless and the fermions are in 
the adjoint representation. The fermionic part of the action is easy to write 
since one has the operator D whose inner fluctuations are 

(10) D A = D + A + JAJ' 1 , A = aj[D, bj] , a h bj GA, A = A* 



using the age of the universe in Planck units to estimate the spatial Euclidean direc- 
tions and the inverse temperature (3 = 1/kT also in Planck units, to set the size of the 
imaginary time component of the Euclidean M. 
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In the Einstein- Yang-Mills system so obtained, all fields involved are mass- 
less. 

We now consider the product M x F of a Riemannian geometry M (of 
dimension 4) by the finite space F of JCO-dimension 6 which was determined 
above. The computation shows that (c/ [9j) 

• The inner fluctuations of the metric give an U(l) x SU(2) x SU(3) 
gauge field and a complex Higgs doublet scalar field. 

• The spectral action Tr(/(D/A)) plus the antisymmetric bilinear form 
(J£, DaX]) on chiral fermions, gives the Standard Model minimally 
coupled to gravity, with the Majorana mass terms and see-saw mech- 
anism. 

• The gauge couplings fulfill the unification constraint, the Yukawa 
couplings fulfill Yi = Ag 2 , where Y2 is defined in eq (jlip . and the 
Higgs quartic coupling also fulfills a unification constraint. 

Most of the new terms occur in the 04 term of the expansion ([8]) . This is the 
case for the minimal coupling of the Higgs field as well as its quartic self- 
interaction. The terms ao and 02 get new contributions from the Majorana 
masses (c/ [9]), but the main new term in ai has the form of a mass term 
for the Higgs field with the coefficient —A 2 . This immediately raises the 
question of the meaning of the specific values of the couplings in the above 
action functional. Unlike the above massless Einstein- Yang-Mills system we 
can no longer take the above action simply as a classical action, would it be 
because of the unification of the three gauge couplings, which does not hold 
at low scale. The basic idea proposed in [8] is to consider the above action 
as an effective action valid at the unification scale A and use the Wilsonian 
approach of integrating the high frequency modes to show that one obtains 
a realistic picture after "running down" from the unification scale to the 
energies at which observations are done. This approach is closely related 
to the approach of Reuter [15], Dou and Percacci [16], [17]. The coarse 
graining uses a much lower scale p which can be understood physically as 
the resolution with which the system is observed. The modes with momenta 
larger than p cannot be directly observed and their effect is averaged out 
by the functional integral. In fact the way the renormalization group is 
computed in [16] shows that the derivative pd p T p of the effective action is 
expressed as a trace of an operator function of the propagators and is thus 
of a similar nature as the spectral action itself, though the trace involves 
all fields and not just the spin ^ fields as in the spectral action. It is an 
open question to compute the renormalization group flow for the spectral 
action in the context of spectral triples. One expects, as explained above, 
that new terms involving traces of functions of the bosonic propagatoijj 
S p SD Tr(f (D / A)) will be generated. The idea of taking the spectral action 
as a boundary condition of the renormalization group at unification scale 



We thank John Iliopoulos for discussions on this point. 
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generates a number of severe tests. The first ones involve the dimensionless 
couplings. These include 

(1) The three gauge couplings 

(2) The Yukawa couplings 

(3) The Higgs quartic coupling 

As is well known, the gauge couplings do not unify in the Standard Model 
but the meeting of g 2 and gs specifies a "unification" scale of ~ 10 17 GeV. 
For the Yukawa couplings the boundary condition gives 

(ii) y 2 = a 9 \ Y 2 = J2(y a v) 2 + (ye) 2 + HyZ) 2 + z(y a d ) 2 - 

17 

This yields a value of the top mass which is 1.04 times the observed value 
when neglecting^ the Yukawa couplings of the bottom quarks etc... and is 
hence compatible with experiment. The Higgs quartic coupling (scattering 
parameter) has the boundary condition of the form: 

A(A) = si ~ dl 

The numerical solution to the RG equations with the boundary value Ao = 
0.356 at A = 10 17 GeV gives \{Mz) ~ 0.241 and a Higgs mass of the order of 
170 GeV. This value now seems to be ruled out experimentally but this might 
simply be a clear indication of the presence of some new physics, instead of 
the "big desert" which is assumed here in the huge range of energies between 
10 2 GeV and 10 17 GeV. To be more precise the above "prediction" of the 
Higgs mass is in perfect agreement with the one of the Standard Model, 
when one assumes the "big desert" (c/ [E]). I n a forthcoming paper [18] 
we show that the choice of the spectral function / could play an important 
role, even when it varies slightly from the cutoff function. This is related 
to the fact that the vev of the Higgs field is proportional to the scale A 
and thus higher order corrections do contribute. This will cause the relation 
between the gauge coupling constants to be modified and to change the 
Higgs potential. Such gravitational corrections are known to cause sizable 
changes to the Higgs mass |20j . 

The next tests involve the dimensionful couplings. These include 

(1) The inverse Newton constant Z g = 1/G. 

(2) The mass term of the Higgs. 

(3) The Majorana mass terms. 

(4) The cosmological constant. 

Since our action functional combines gravity and the Standard Model, the 
analysis of [16] applies, and the running of the couplings Z which have 
the physical dimension of the square of a mass is well approximated by 
f3 z = a\k 2 where the parameter k is fixing the cutoff scale but is considered 
itself as one of the couplings, while the coefficient a± is a dimensionless 



'See [9] for the precise satement. 
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number of order one. For the inverse Z g of the Newton constant, one gets 
the solution: 

(12) Z g = Z g (l + \a 1 1 t) 

which behaves like a constant and shows that the change in Z g is moderate 
between the low energy value Z g at k = and its value at k = mp the Planck 
scale, for which -|- = 1. We have shown in [9] that a relation between the 
moments of the cutoff function / involved in the spectral action, of the form 
/2 ~ 5/o suffices to give a realistic value of the Newton constant, provided 
one applies the spectral action at the unification scale A ~ 10 17 GeV. The 
above discussion of the running of Z g shows that this yields a reasonable 
low energy value of the Newton constant G. 

The form 0z = a\k 2 of the running of a coupling with mass 2 dimension 
implies that, as a rule, even if this coupling happens to be small at low 
scale, it will necessarily be of the order of A 2 at unification scale. For the 
Majorana mass terms, we explained in [9] why they are of the order of 
A 2 at unification and their role in the see-saw mechanism shows that one 
should not expect them to be small at small scale, thus a running like (fT2j) 
is realistic. Things are quite different for the mass term of the Higgs. The 
spectral action delivers a huge mass term of the form —A 2 H 2 and one can 
check that it is consistent with the sign and order of magnitude of the qua- 
dratic divergence of the self-energy of this scalar field. However though this 
shows compatibility with a small low energy value it does by no means al- 
low one to justify such a small value. Giving the term —A 2 if 2 at unification 
scale and hoping to get a small value when running the theory down to low 
energies by applying the renormalization group, one is facing a huge fine 
tuning problem. Thus one should rather try to find a physical principle to 
explain why one obtains such a small value at low scale. In the noncommu- 
tative geometry model M x F of space-time the size of the finite space F 
is governed by the inverse of the Higgs mass. Thus the above problem has 
a simple geometric interpretation: Why is the space F so largefl in Planck 
units? There is a striking similarity between this problem and the problem 
of the large size of space in Planck units. This suggests that it would be 
very worthwhile to develop cosmology in the context of the noncommutative 
geometry model of space-time, with in particular the preliminary step of the 
Lorentzian formulation of the spectral action. 

This also brings us to the important role played by the dilaton field which 
determines the scale A in the theory. The spectral action is taken to be a 
function of the twisted Dirac operator so that D 2 is replaced with e~^D 2 e~^ . 
In [10] we have shown that the spectral action is scale invariant, except for 
the dilaton kinetic energy. Moreover, one can show that after rescaling the 



^by a factor of 10 
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physical fields, the scalar potential of the theory will be independent of the 
dilaton at the classical level. At the quantum level, the dilaton acquires 
a Coleman- Weinberg potential [21] and will have a vev of the order of the 
Planck mass |22j. The fact that the Higgs mass is damped by a factor of 
e~ 2< ^, can be the basis of an explanation of the hierarchy problem. 
In this paper we investigate the accuracy of the approximation of the spec- 
tral action by the first terms of its asymptotic expansion. We consider the 
concrete example given by the four-dimensional geometry x Sp where S% 
is the round sphere of radius a as a model of space, while Sp is a circle of ra- 
dius j3 viewed as a model of imaginary periodic time at inverse temperature 
p. We compute directly the spectral action and compare it with the sum of 
the first terms of the asymptotic expansion. In section two we start with the 
round sphere and use the known spectrum of the Dirac operator together 
with the Poisson summation formula, to estimate the remainder when using 
a smooth test function. This is then applied to the four-dimensional space 
&a x where it is shown that, for natural test functions, the spectral action 
is completely determined by the first two terms, with an error of the order 
of 10 -<J where a is the inner diameter A//, fi = inf (a, p) in units of the 
cutoff A. Thus for instance an inner diameter of 10 in cutoff units yields the 
accuracy of the first hundred decimal places, while an inner diameter of 10 31 
corresponding to the visible universe at inverse temperature of 3 Kelvin and 
a cutoff at Planck scal^l, yields an astronomical precision of 10 62 accurate 
decimal places. This is then extended in the presence of Higgs fields. The 
above direct computation allows one to double check coefficients in the spec- 
tral action. It also implies, for x S\, the vanishing of all the Seeley-De 
Witt coefficients a,2n, n > 2, in the heat expansion of the square of the Dirac 
operator. This is confirmed in section three, by a local computation of the 
heat kernel expansion, where it is shown that 04 and vanish due to subtle 
cancelations. 



2. Estimate of the asymptotics 

The number A" (A) of eigenvalues of \D\ which are < A 

(13) AT (A) = # eigenvalues of D in [-A, A], 

is a step function A" (A) which jumps by the integer multiplicity of an eigen- 
value whenever A belongs to the spectrum of \D\. This integer valued func- 
tion is the superposition of two terms, 

N(A) = (AT(A))+AT osc (A). 

The oscillatory part AT osc (A) is generically the same as for a random matrix. 
The average part (A"(A)) is computed by a semiclassical approximation from 



while the age of the universe in Planck units gives Aa ~ 10' 
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local expressions involving the familiar heat equation expansion and will now 



be carefully defined assuming an expansion of the for 

(14) Trace (e" iA ) ~ ^ a a t a (t -» 0) 
for the positive operator A = D 2 . One has, 

(15) A~ s / 2 = — ^r- / e-^t^-Ut 

and the relation between the asymptotic expansion (|14p and the £ function, 

(16) Cd(s) = Trace (A" s / 2 ) 

is given by, 

• a < gives a po/e at —2a for £d with 

(17) Res s= _2aCr>(s) = tJ— 

r(-or) 

• a = (no logi term) gives regularity at for £d with 

(18) ( D (0) = ao . 

For simple superpositions of exponentials, as Laplace transforms, 



CO 



(19) /(«) = / e- sn /i(s)ds 
we can write formally, 

/>oo 

(20) f(tA) = / e~ s * A h(s) ds 

Jo 

and 

/■oo 

(21) Trace (/(tA)) ~ V a a i a / s a /»(s) (is . 

./o 

For a < one has, 

1 /'OO 



T(-a) 
and 



s «/ l ( a )d s = _ -/ f{v)v- a - l dv 

r(-a) 7o 



so that 



1 f 00 
Trace (/(tA)) ~ V - Res s= _ 2a Q D {s) / /(«) t^ 1 du t c 

(22) + Cd(0) /(0) + Y j a a t a s a h(s) ds . 



^the a a defined here is equal to the Seeley-de Witt coefficients a n+ 2 a in dimension 
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Now we assume that the only a > for which a a ^ are integers and note 
that 

/•oo 

(23) / s n h(s)ds = (-lTf^(0), 



so that all the terms a a for a > have vanishing coefficients when / is 
a cutoff function which is constant equal to 1 in a neighborhood of 0. To 
define the average part we consider the limit case f(v) = 1 for < 1 and 
elsewhere and get for the coefficients of ([221 

i poo j.a 

(24) - / f(v)v~ a - l dvt a - 



2 Jo W (-2a) ' 

which, with t = A -2 , gives the following definition for the average part 

(25) (N(A)) -=J2~k ReSs=k Cd{s) + Cd(0) • 

fc>0 

To get familiar with this definition we shall work out its meaning in a simple 
case, 

Proposition 1. Assume that Spec DcZ and that the total multiplicity of 
{±n} is P (n) for a polynomial P(x) = Yl Ck x ■ Then one has 

(N(A))= P(u)du + c, c = y^c k C{-k), 
Jo 

where £ is the Riemann zeta function. 

Proof. One has by construction, with P(x) = ^ CkX k , 

Cd(s) = ]T P(n) n~ s = ^ c fc C(« " k) 

Thus 

Res s=fc Ci)(s) = c fe _i 



and 

k 



fc>0 



The constant Cd(O) is given by 

and is independent of A. □ 
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2.1. The sphere S 4 . We check the hypothesis of Proposition [T] for a round 
even sphere. We recall (|26j) that the spectrum of the Dirac operator for 
the round sphere S n of unit radius is given by 

71 

(26) Spec(.D) = {±(- + k) | k € Z,k > 0} 

where the multiplicity of + k) is equal to 2W ( fc+ /? _1 ) • Thus for n = 4 one 

gets that the spectrum consists of the relative integers, except for { — 1, 0, 1}. 

The multiplicity of the eigenvalue m is 4( fc ^ 3 ) for k + 2 = m which gives, 

for the total multiplicity of ±m 

4 4 3 

P(m) = ~{m + l)m(m — 1) = — (m — m) 

which shows that one gets the correct minus sign for the scalar curvature 
term after integration using Proposition [TJ Thus one gets (up to the nor- 
malization factor | ) 

(27) Tr(| J Dr s ) = C(s-3)-C( S -l) 

This function has a value at s = given by 

x ,, x 1 1 11 
C(-3 - C(-i) = — + — = — 
sv ; sv ; 120 12 120 

which, taking into account the factor | from normalization, matches the 
coefficient x 4 which appears in the spectral action in front of the Gauss- 
Bonnet term, as will be shown in §3) 

2.2. The sphere S 3 . We now want to look at the case of S 3 and determine 
how good the approximation of ()22[) is for test functions. 

In order to estimate the remainder of (1221) in this special case we shall use 
the Poisson summation formula 

(28) ^M^) = E^ n )' ^0= / h ( l 

z z R 

or rather, since the spectrum is | + Z in the odd case, the variant 

(29) J>(n+^) = J>imn) 

z z 

(obtained from (f28j) using fo(u) = + ^)). 

In the case of the three sphere, the eigenvalues are ±(| + k), for k > with 
the multiplicity 2( fe ^ 2 ). Thus re + i has multiplicity n(n + 1). This holds 
not only for ra > but also for n € Z since the multiplicity of — (n + |) 
is n(n + 1) = m(m + 1) for m = — n — 1. In particular ±^ is not in the 
spectrum. Thus when we evaluate Tr(/(D/A)), with / an even function, we 
get the following sum 

(30) Tr(/(.D/A)) = £ n(n + l)/((n + \)/K) 



<u)e- 2mxu du 
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We apply (|29j) with g(u) = (u 2 — \)f{u/A). The Fourier transform of g is 

g(u)e- 2nixu du = [ {u 2 - -)f{u/K)e- 2 ™ xu du 
Jr 4 

= A 3 f v 2 f(v)e- 2mAxv dv - ]A f f(v)e- 2mAxv dv 
Jr 4 J R 

We introduce the function which is the Fourier transform of v 2 f(v) and 
we thus get from ([29]) . 

(31) Tr(f(D/A)) = A 3 ^(-l)"/( 2 )(An) - ±A£(-l) B /(An) 

z z 

If we take the function / in the Schwartz space <S(R), then both / and 
have rapid decay and we can estimate the sums 

£ |/(An)| < C k A- k , l/ (2) (An)| < C k A- k 

which gives, for any given k, an estimate for a sphere of radius a of the form: 

(32) Tr(/(D/A)) = (Aa) 3 / v 2 f(v)dv - - (Aa) / f(v)dv + 0((Aa)- fc ) 

Jr 4 Jr 

The radius simply rescales D and enters in such a way as to make the product 
Aa dimensionless. This can be seen by noting that the ratio contains the 
term TeSf^ and the radius enters as ^ in the inverse dreibein e«. Note 
that, provided that > 1 one controls the constant in front of (Aa)~ k from 
the constants Cj with 

\x k f(x)\ <c u \x k f 2 \x)\<c 2 . 

To get an estimate of these constants Cj, say for k = 2, one can use the 
L 1 -norms of the functions Af(v) and A(v 2 f(v)) where A = — d 2 is the 
Laplacian. If we take for / a smooth cutoff function we thus get that the Cj 
are of order one. 

In fact we shall soon get a much better estimate (Corollary U] below) which 
will show that, for suitable test functions, a size of N in cutoff units, Aa ~ N, 
already ensures a precision of the order of e~ N2 . We shall work directly with 
the physically more relevant model consisting of the product S s x S 1 viewed 
as a model of the imaginary time periodic compactification of space-time at 
a given temperature. Our estimates will work well for a size in cutoff units 
as small as ~ 10 and will give the result with an astronomical precision 
for larger values. These correspond to later times since both the radius of 
space and the inverse temperature are increasing functions of time in this 
simple model. 
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2.3. The product S 3 x S 1 . We now want to move to the 4-dimensional 
Euclidean case obtained by taking the product M = S 3 x S 1 of S 3 by a 
small circle. We take the product geometry of a three dimensional geometry 
with Dirac operator D% by the one dimensional circle geometry with Dirac 

(33) Di=iiV fl 
so that the spectrum of D\ is jj(Z + |). 

Lemma 2. Lei D 6e i/ie Dirac operator of the product geometry 

m N n _( L> 3 ® 1 + * ® £>i 

1 j y £> 3 <g> 1 - % ® I>i 

T/ie asymptotic expansion for A — > oo o/ i/ie spectral action of D is given by 

(35) Tr(h(D 2 /A 2 )) ~ 2/?ATr(A:( J D|/A 2 )), 
where the function k is given by 

poo 

(36) fc(a;) = / (it - x)~ 1/2 



Proof. By linearity of both sides in the function h (using the linearity of 
the transformation (|36p ) it is enough to prove the result for the function 
h(x) = e~ bx . One has 



D 2 



Dl ® 1 + 1 (8) £> 2 

£>§ (8) 1 + 1 (8) £>? 



and 

Tr( e - bD2 / A2 ) = 2Tr(e- feD ?/ A2 )Tr(e- feD a 2 / A2 ) 

Moreover by ([33]) the spectrum of Z>i is -^(Z + |) so that, using ([29]) . and 
for fixed (3 and 6, one has for all k > 0, 

Thus 

Tr(e _ bD 2 /A2) = 2/3ATr( ^ Fr l/2 e -6Cl/A^ ) + 0(A -fc+3 ) 
and the equality (|35p follows from 

/>oo 

/ (u - x)- 1 ' 2 e~ bu du = ^b- l > 2 e~ bx 

J X 

which shows that the function k associated to h(x) = e~ bx by the linear 
transformation (|36p is k(x) = ^/tt b~ 1 / 2 e~ bx . □ 



One can write (J36J) in the form 

POO 

(37) k{x) = / v~ l/2 h(x + v)dv, 
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which shows that k has right support contained in the right support of h 
i.e. that if h vanishes identically on [0,00 [ so does k. It also gives a good 
estimate of the derivatives of k since 

POO 

d%k(x)= / v- l/2 d™h{x + v)dv. 



In fact, in order to estimate the size of the remainder in the asymptotic 
expansion of the spectral action for the product M = S 3 x S 1 , we shall now 
use the two dimensional form of (1291). 



(38) £ g(n + \ , m + \) = £(-l)«+^(n, m) 

T 2 Tt? 

where the Fourier transform is given by 

(39) 9{x,y)= [ g{u,v)e~ 2m{xu+yv Uudv 



For the operator D of (|34|) . and taking for -D3 the Dirac operator of the 
3-sphere S 3 of radius a, the eigenvalues of D 2 / A 2 are obtained by collecting 
the following 

(i + n) 2 (Aa)" 2 + (^ + m) 2 (A/3)- 2 , n,mEZ 
with the multiplicity 2n(n + 1) for each n, m E Z. Thus, more precisely 

Tr{h(D 2 /A 2 )) = Mn + + ™) 2 (Aa)~ 2 + d + ™? W % ) 

1? 

which is of the form: 

(40) Tr(h(D 2 /A 2 )) = J2 9(n + ~ m + ~) 

1? 

where 

(41) g(u, v) = 2{u 2 - ^)h{u 2 (Aa)~ 2 + v 2 (A/3)^ 2 ) 
One has 

£(0, 0) = / g(u, v)dudv = 2 [ (u 2 - \)h(u 2 {Aa)~ 2 + v 2 (A(3)~ 2 )dudv 
Jr 2 Jr 2 4 

= 2(Aa)(A/3) / ((Aa) 2 x 2 - -)h(x 2 + y 2 )dxdy 
Jr 2 4 

using u = x (Aa) and v = y (A/3) . Thus we get: 

roo poo 

(42) 3(0, 0) = 2ir (A/3) (Aa) 3 / h(p 2 )p 3 dp - vr (A/3) (Aa) / h(p 2 )pdp 

Jo Jo 

To estimate the remainder, given by the sum 

(-l) n+m g(n,m) 

(n,m)^(0,0) 
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we treat separately the Fourier transforms of 

gi{u,v) = u 2 h(u 2 {Aa)- 2 + v 2 ((3A)~ 2 ) , ga(u,«) = h(u 2 {A a y 2 + v 2 (A/3)" 2 ) 
One has 

g 2 (n,m) = [ g 2 {u,v)e- 2 ^ nu+mv Uudv 
Jr 2 

= A 2 f3a [ h(x 2 + y 2 )e- 2m( - nAax+mA ^dxdy = A 2 (3an 2 (nAa, mA/3) 
Jr 2 

where the function of two variables k 2 (u,v) is the Fourier transform, 

(43) k 2 (u,v)= f h(x 2 + y 2 )e- 2 ^ ux+vy) dxdy = k(u 2 + v 2 ) 

Jr 2 

The function k is related to the function k(x) defined by (|36p . and one has 

(44) k{u 2 ) = ( k(x 2 )e- 27Tiux dx 

Jr 

so that k(u 2 ) is the Fourier transform of k(x 2 ). 
For gi one has, similarly, 

gi (n,m)= [ gi { Ul v)e- 2ni(nu+mv) dudv 
Jr 2 

= A 4 /3a 3 ! x 2 h(x 2 + y 2 )e- 2ni{nAax+mAI3y) dxdy = A 4 /3a 3 Kl (nAa, mA/3) 
Jr 2 

where the function of two variables k,\ (u, v) is the Fourier transform, 

ki («,«)= / x 2 h(x 2 + y 2 )e- 2ni{ux+vy) dxdy 
Jr 2 

which is given in terms of (|43[) by 

(45) v) = -7t- 2 (u 2 k"(u 2 + v 2 ) + \'{u 2 + v 2 )) 

Now for any test function h in the Schwartz space <S(R), the function 
x 2 h(x 2 + y 2 ) is in the Schwartz space 5(M 2 ) and thus we have for its Fourier 
transform, and any k > 0, an estimate of the form 

(46) |«i(«,«)| <C k (u 2 + v 2 )- k 
We thus get, for k > 2, 

I £ (-ir +m 5i(™,m)| < £ l5i(n,m)| 

(n,m)^(0,0) (n,m)^(0,0) 

= A 4 /3a 3 £ |Ki(raAa,mA/?)| < C fe A 4 /?a 3 £ ((nAa) 2 +(mA/?) 2 )- fc 

(ra,m)^(0,0) (n,m)^(0,0) 

< C fe A 4 /?a 3 (A / u)- 2fc £ (n 2 + m 2 )" fc , /x = inf (a, /5) 

(n,m)^(0,0) 

We thus get, using a similar estimate for g 2 , 
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Theorem 3. Consider the product geometry x S\. Then one has, for 
any test function h in the Schwartz space <S(K), the equality 

roo roc 

(47) Tr(h(D 2 /A 2 )) = 2vrA 4 /?a 3 / h(p 2 )p 3 dp-7rA 2 (3a / h{p 2 )pdp+e{A) 

Jo Jo 

where e(A) = 0(A~ k ) for any k is majorized by 

|e(A)| < 2A 4 /3a 3 ^ l«i(™Ao, mA/3)|+^A 2 /3a ^ |/e 2 (nAa, mA0)\ . 

(n,m)^(0,0) " (n,m)^(O,0) 

with Kj defined in (|43f) and (|45p . 

This implies that all the Seeley coefficients a2n vanish for n > 2, and we 
shall check this directly for 04 and a§ in £j3j 

This vanishing of the Seeley coefficients does not hold for the 4 sphere and it 
is worth understanding why one cannot expect to use the Poisson summation 
in the same way for the 4 sphere. The problem when one tries to use the 
Poisson formula as above is that, e.g. for the heat kernel, one is dealing with 
a function like |x|e -fcr which is not smooth and whose Fourier transform 
does not have rapid decay at 00. 

2.4. Specific test functions. We shall now concretely evaluate the re- 
mainder in Theorem [3] for analytic test functions of the form 

(48) h{x) = P{vx)&-™ 

where P is a polynomial of degree d. The Fourier transforms of the functions 
of two variables x 2 h(x 2 + y 2 ) and h(x 2 + y 2 ) are of the form 

Kj (u,v) = Pj{u,v)e-^ u2+v ^ 
where the Pj are polynomials. More precisely, since the Fourier transform of 

\ I 2 1 2\ 1 (u 2 +v 2 ) 

e v*^ +v ) is je~ n * one obtains the formula for P2 by differentiation 
at A = 1 and get 

1 _ (u 2 +v 2 ) 

k 2 (u,v) = P(-d x )\=ije n a 

which is of the form 

«2(«, v) = Q(tt(u 2 + v 2 ))e-^ u2+v ^ 

where Q is a polynomial of degree d. The transformation P 1— > Q = T(P) is 
given by 

(49) Q(z) = P(-d x ) x=1 ~e-i 
Moreover one then gets 

ki(u,v) = -(2ir)~ 2 dlK 2 (u,v) 

= (u 2 Z 1 (tt(u 2 + v 2 )) + Z 2 (tt(u 2 + v 2 )))e-^ u2+v ^ 
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where 

(50) Z 1 = -Q + 2Q' -Q" , Z 2 = ^-(Q-Q') 

We let Cp be the sum of the absolute values of the coefficients of Q = T(P). 



Corollary 4. Consider the product geometry S 3 x S\. Let p = inf(a, j3). 
Then one has, with h any test function of the form (|48p , the equality 

poo poo 

(51) Tr(/i(D 2 /A 2 )) = 2vrA 4 /3a 3 / h{p 2 )p 3 dp-TrA 2 f3a / h(p 2 )pdp+e(A) 

Jo Jo 

where, assuming pA > \/d(l + logd) and pA > 1, 

(52) |e(A)| < Ce-i^ 2 , C = A 4 /3a 3 C P (8 + 6d + 2d 2 ) 
Proof. One has 

x k e -x/2 < x ^ Vx > 3^ + log fc ) 

Thus, for (n, m) 7^ (0, 0) one has 

\K 2 {nAa,mAp)\ < c pe -^ nAa ^ +{mA ^ 

since 7r((nAa) 2 + (mA/5) 2 ) > 3d(l + logd). Moreover, since e~5^ A ) 2 < i, 
one gets 

£ e -|((nAa) 2 +(mA/3) 2 ) < gg-f (M) 2 
(n,m)^(0,0) 

and 

I k 2 (nAa, mA/3) I < 8 C P e" i (M)2 

(ra,m)^(0,0) 

A similar estimate using (I50D yields 

\Ki{nAa,mA[3)\ < (2 + 3d + d 2 ) Cpe~i (M)2 

(n,m)^(0,0) 

Thus by Theorem El the inequality ([52]) holds for 

C = C P (2A 4 /3a 3 (2 + 3d + d 2 ) + 4A 2 /3a) . 
One then uses the hypothesis pA > 1 to simplify C. □ 

The meaning of Corollary U] is that the accuracy of the asymptotic expansion 
is at least of the order of e _ "2^ A - ) . Indeed the term A 4 /?a 3 in (f52|) is the 
dominant volume term in the spectral action and the other terms in the 
formula for C are of order one. Thus for instance for a size pA ~ 100 one 
gets that the asymptotic expansion accurately delivers the first 6820 decimal 
places of the spectral action. Note that some test functions of the form (I48p 
give excellent approximations to cutoff functions, in particular 

(53) MaO = En?e-™ 
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The graph of h n (x 2 ) is shown in Figure [T] for n = 20. For h = /2.20 the 
computation gives Cp(8 + 6d + 2d 2 ) < 2 x 10 6 so that this constant only 
interferes with the last six decimal places in the above accuracy. 

In our simplified physical model we test the approximation of the spectral 
action by its asymptotic expansion for the Euclidean model 

where space at a given time t is given by a sphere with radius a(t) and (5{t) 
is a uniform value of inverse temperature. One can then easily see that the 
above approximation to the spectral action is fantastically accurate, going 
backwards in time all the way up to one order lower than the Planck energy. 
In doing so the radius a(t) varies between at least ~ 10 61 Planck units and 
10 Planck units (i.e. 10~ m), while the temperature varies between 2.7° K 
and (l0 31 )°-RT. It is for an inner size less than 10 in Planck units that the 
approximation does break down. 

Remark 5. For later purpose it is important to estimate the constant Cp 
in terms of the coefficients of the polynomial P. Let then P{z) = z n . One 
has h(x) = (■Kx) n e~ 7rx and the function k(x) associated to h by (f36l) is 

k(x)= [ h{x + y 2 )dy = Ti n e- 7Tx y" ("' S ]x n - k [ y 2k e- ny2 dy 
Jr W Jr 

^ ' 

To obtain Q = T(P) one then needs to compute the Fourier transform n(u 2 ) 
of the function k(x 2 ) as in (|44p . The Fourier transform of (Trx 2 ) m e~' KX is 

£ m (u) = (-4Tr)- m d 2 u m e-™ 2 = L m (v™ 2 )e-™ 2 

and one checks, using the relation 

L m+1 (z) = 1/2((1 - 2z)L m {z) + (-1 + Az)L' m {z) - 2zL" m {z) 

that the sign of the coefficient of z k in L m (z) is (— l) fc . Thus the sum of the 
absolute values of the coefficients of L m is equal to L m (— 1) = £ m (i7r -1 / 2 )/e. 
Thus since the above sum giving k(x) has positive coefficients we get that, 
for P(z) = z n , the constant Cp is given by Q(tt(u 2 + v 2 )e~ w ( u +v ))/e for 
(u, v) = (iTr" 1 / 2 , 0) , which gives 

C P = [ ir n (y 2 + x 2 ) n e- wy2 -™ 2+2 ^ x ~ 1 dxdy. 

One then gets 

/•oo 

(54) C P <2 u 2n+1 e-( u ~V 2 du = 0{\ n n\) , VA > 1 . 

Jo 
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id 



Thus, for an arbitrary polynomial P(z) = Y2o a k z ^ one nas 

poo 

(55) C P <2 \P\(u 2 )e- {u - 1)2 udu, \P\(z) =J2\a k 

Jo 



2.5. The Higgs potential. We now look at what happens if one performs 
the following replacement on the operator 

D 2 ^ D 2 + H 2 

where H is a constant. This amounts with the above notations to the 
replacement 

(56) h{u) >-> h(u + H 2 /A 2 ) . 

As long as H 2 /A 2 is of order one, we can trust the asymptotic expansion 
and we just need to understand the effect of this shift on the two terms of 
(|47p . We look at the first contribution, i.e. 

poo poo 

27rA 4 /3a 3 / h(p 2 )p 3 dp = ^A 4 /3a 3 / uh(u)du 
Jo Jo 

We let x = H 2 /A 2 , and get, after the above replacement, 

' poo poo px 

uh(u + x)du = I (v — x)h(v)dv = / (y — x)h{v)dv— I {y — x)h{v)dv 
Jx Jo Jo 

poo poo px 

vh(v)dv — x h{v)dv — I (v — x)h(v)dv 
io Jo Jo 

The first term corresponds to the initial contribution of 7rA 4 /3a 3 uh(u)du. 

The second term gives 

poo poo 

(57) - TiA 4 (3a 3 x / h(v)dv = -irA 2 f3a 3 H 2 / h(v)dv 

Jo Jo 

which is the expected Higgs mass term from the Seeley-de Witt coefficient 
a>2. To understand the last term we assume that h is a cutoff function. 

Lemma 6. // h is a smooth function constant on the interval [0, c], then 
for x = H 2 / A 2 < c the new terms arising from the replacement (|56p are 
given by 

r-oo i ^ 

(58) - 7rA 2 /3a 3 / h{v)dv H 2 + -ir[3ah(0) H 2 + -irpa 3 h(0) H 4 

Jo 2 2 

Proof. For the perturbation of it A 4 (3a 3 / °° uh(u)du, besides ([57]) . we just 
need to compute the last term — Jq(v — x)h{v)dv, and one has 

px px 

- (v- x)h(v)dv = h(0) (x- v)dv = -h(0)x 2 
Jo Jo 2 

since h is constant on the interval [0,x]. 
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We then look at the effect on the second contribution, i.e. 

roc j roc 

-TrA 2 /3a h(p 2 )pdp = --vrA 2 (3a h(u)du 
Jo * Jo 

We let, as above, x = H 2 /A 2 , and get 

roc roc rx 

h(u + x)du = / h[v)dv = I h(v)dv — / h(v)dv 



Thus the perturbation, under the hypothesis of Lemma [6] is 
-^7rA 2 f3a(-xh(0)) = -7r/3a/i(0) H 2 

□ 

The three terms in formula (I58p correspond to the following new terms for 
the spectral action 

• The Higgs mass term coming from the Seeley-de Witt coefficient 02- 

• The RH 2 term coming from the Seeley-de Witt coefficient 04. 

• The Higgs potential term in P 4 coming from the Seeley-de Witt 
coefficient 04. 

We can now state the analogue of Theorem [3] as follows 

Theorem 7. Consider the product geometry S% x S\. Let p = inf(a,/3). 
Then one has, with h any test function of the form (1481) , the equality 



roc roc 

Tr(h((D 2 + H 2 )/A 2 )) = 2vrA 4 /3a 3 / h(p 2 )p 3 dp - nA 2 (3a / h(p 2 )pdp 

Jo Jo 

+irA 4 pa 3 V(H 2 /A 2 ) + \vA 2 f3a W{H 2 /A 2 ) + e(A) 

where 

roc rx 

(59) V{x)= I u{h{u + x) - h(u))du , W(x) = / h{u)du 

Jo Jo 

and, assuming pA > \/ d(l + log d), pA > 1, and H 2 A~ 2 < c/tt, 

(60) |e(A)| < Ce'^V 2 , C = A 4 /3a 3 C^(8 + 6d + 2d 2 ) 
where, with P(z) = Ylo a k zk ° ne h as 

roc 

C' P = A \P\{u 2 + c)e-( u -V 2 udu, \P\(z) = y2\a k \z k 
Jo 



Proof. The new terms simply express the replacement (I56P in the formula of 
Theorem [3l The new function h thus obtained is still of the form (|48p since 
it is obtained from h by a translation. It thus only remains to estimate Cp 
where P is the polynomial such that h(u) = P(7ru)e _7r " . For P(z) = z n the 
constant Cp for a translation u ^ u + x, x>0of the variable, is less than 
the constant Cp x for the polynomial 



irx) n - k Uu) k 
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Thus, by Remark [5l (i55j) . the constant Cp x is estimated by 



/•oo 

C Px <2 (u 2 + ■Kx) n e- {u ~ l) \du 
Jo 



which is an increasing function of x and thus only needs to be controlled for 
x = c/tt in our case. □ 

For instance, for h = /120 the computation gives C' P (8 + 6d + 2d 2 ) < 3 x 10 7 
for c = 1, so that this constant only interferes with the last seven decimal 
places in the accuracy which is the same as in Corollary 01 
Moreover as shown in Lemma [6l when h is close to a true cutoff function 

(61) vrA 4 /5a 3 V(H 2 /A 2 ) + \^k 2 (5a W(H 2 /A 2 ) 

f°° 1 1 

= -2nA 2 Pa 3 / h(p 2 )pdpH 2 + -vr/3a/i(0) H 2 + -irpa 3 h(0) H 4 + 6 
Jo 2 2 

where the remainder 5 is estimated from the Taylor expansion of h at 0. For 

instance for the functions h n of (|53p . one has by construction < h n (x) < 1 

for all x and since 

h n (x) = 1 - k)x n+k+1 , a(n, k) = (-l) k /({n + k + l)n\kl) 

one gets, for h = h n the estimate 

n+3 rj.n+2 

l J l ^ ttAW , f + ^A 2 /3a-^— - , x = i? 2 /A 2 . 

(n + 3)(n + lj! 2(n + 2)! 

While the function W is by construction the primitive of h, and is increasing 
for h > one has, under the hypothesis of positivity of h, 

Lemma 8. The function V{x) is decreasing with derivative given by 

roo 

V'(x) = - h(v)dv 

J x 

The second derivative ofV(x) is equal to h(x). 
Proof. One has 

/'OO roo 

V'(x)= / uh'(u + x)du = [uh(u + x)]^ - / h(u + x)du 
Jo Jo 

which gives the required results. □ 

3. Seeley-De Witt coefficients and Spectral Action on S 3 x S 1 

In this section we shall compute the asymptotic expansion of the spectral 
action on the background geometry of S 3 x S 1 using heat kernel methods. 
This will enable us to check independently the accuracy of the estimates 
derived in the last section. This background is physically relevant since it 
can be connected with simple cosmological models. We refer to [23], [8] for 
the formulas and the method of the computation. The general method we 
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use is also explained in great detail in a forthcoming paper [To] . We start 
by computing ao : 

s 3 s 1 



' (2vr 2 a 3 ) (2vr/3) = tt/3. 



4vr 2 

where (5 is the radius of Sp and the volume of the three sphere of radius 
a is 2vr 2 a 3 [25] . 
Next we calculate 02 

where E is defined from the relation 



D 2 = — {gTVuVu + E) 



where for pure gravity we have 



so that (using Tr(l) =4 



E = -\r 
4 



02 = h ("if) 

since the curvature is constant. The curvature tensoi[3 is, using the coordi- 
nates of [25] for the three sphere S% with labels i, j, k, I and the label 4 for 
the coordinate in Si, 

Rijki = ~a~ 2 {gikQji ~ 9il9jk) , i,3, k, I = 1, • ■ • 3 
Rijki = 
RiijA = 

where gij is the metric on the three sphere as in [25J. The Ricci tensor is 
given, following the sign convention of [23] which introduces a minus sign in 
passing from the curvature tensor to the Ricci tensor, by 

Rij = —g kl Rikji = 2a 2 gij 
Ra = 

i?44 = 

Thus the scalar curvature is 

6 



R = g lJ R. 



12 the sign convention for this tensor is the same as in [23] 
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and the a 2 term in the heat expansion simplifies to 

T 



a 2 = -7r (3 'a , o 



Next for 04 we have 



M 

+ 2R llvpa R^ upa + 60i2E + 180£ 2 + 60£ ;M " + 30 ft^O^) 
where for the pure gravitational theory, we have 

E = _ 4^' Qfj,v = -^R^ u ah lab 
In this case it was shown in [8] that 04 reduces to 



2 1 



° 4 = i 36o/ ^ " " ^ 
(62) = / ^^(-lSC 2 ^ + 

which is obviously scale invariant. The Weyl tensor Cavpa is defined by 

Cfjippa — Rpvpcr ~\~ — (Rpp9va RupdpXT Rpa9vp ~\~ Rvo9 'pp) 
- g {9pp3va - 9vp9p(j) R 

This tensor vanishes on S 3 x S 1 as can be seen by evaluating the components 

Cijki = a' 2 [-(gik9jl ~ giWjk) + ^{9ik9jl ~ 9il9jk) ~ (9ik9jl ~ 9il9jk)] = 
Cijki = 

Cj4fc4 = 

Similarly the Gauss-Bonnet term 

R*R* = -e^po ' eafrsR^ 011 * R pa l5 

= 

The next step of calculating a6 is in general extremely complicated, but 
for spaces of constant curvature the expression simplifies as all covariant 
derivatives of the curvature tensor, Riemann tensor and scalar curvature 
vanish. The non-vanishing terms are, using Theorem 4.8.16 of [23] and the 
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above sign convention for the Ricci tensor Rp V and the scalar curvature, 

a6 = ihl dAx ^ Ti (9^7! ( 35jRS _ 42Rr2 ^ + A2RR l»p° 

208Rf lu Rf 1 pR U p \Q1RppR vr jRp V p a ^iRppRppanRvpan 

4:ARp u p (T Rp l , K xRp aK \ 80Rp,L>paRpKp\RvK(T\) 
+ ggQ ( ^^^HU^Up^pfl ^RpVpCT^pU^pcJ ^RpU^pp^-Up + ^R^pu 

+60E 3 + 30EnJ„ + 30RE 2 + 5R 2 E - 2R\ V E + IR^^E)) 

We can now compute each of the above eighteen terms. These are listed in 
an appendix. Collecting these terms we obtain that the integrand is 

4 a -6 

- -^-^ (-35 • 6 3 + 42 • 72 - 42 • 72 + 208 • 24 - 192 • 24 + 48 • 24 - 44 • 24 - 80 • 6) 

4a- 6 f 15-27 5-27 

9 + 18 - 12 + 45 + — 15-27 + 10-27- 36 + 36 



360 



= •-*(§-§) = " 

implying that 

a 6 = 0, 

which shows that the cancelation is highly non-trivial. We conclude that 
the spectral action, up to terms of order is given by 

/*oo /*oo / 1 \ 

S = A 4 J xh (x) dx (irpa 3 ) - A 2 J h (x) dx 1 ir/3a-\ + O (A" 4 ) 



After making the change of variables x = p we get 



S = (tt/3A) 



2 (Aa) 3 / p 3 h (p 2 ) dp - (Aa) / ph (p 2 ) dp 



+ 0(A" 



This confirms equation (|47p and shows that, to a very high degree of accu- 
racy, the spectral action on S 3 x 5 1 is given by the first two terms. 

Remark 9. It is worth noting that one can also check the value of the 
Gauss-Bonnet term on S and show that it agrees with the value obtained 
in (|27p . To see this note that the Riemann tensor in this case is given by 

M) 

Rp,i>pa = O- {9p,p9u<T 9p,(r9vp) 

which implied that 

Cpi/pa — 
R*R* = 6a~ 4 



^One can double check the value of R*R* using the Gauss-Bonnet Theorem. 
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and thus 



The volume of 5 4 is 



Va 



and this implies that 



1 11 
4^ 60~' 



d A x^/g 



,S' 4 



2vri 4 8^ 2 4 

-a = a 

3 



r(|) 



ii , 

04 = X 4 

360 



which agrees exactly with the calculation of ([27]) based on zeta functions. 

Appendix 

In this appendix we compute the eighteen non- vanishing terms that appear 
in the a§ term of the heat kernel expansion. Using the properties 



R 



Ri = 12a 



35ir 



-A2RR% 

—208Rp U Rp P R up 
-\92RppR UIJ Rp l/ p a 



A8Rp, u RppaK RypaK 



35(6) 3 cT 6 
-42(6) (12) a" 6 

-208 (2) 3 gij g ik g jk = 208 (2) 3 (3) a" 6 
— \92RikRjiRijki 

192 (2) 2 ftfc^ (ftfcOj; - fi^jfe) oT 6 

192 (24) a -6 

—48 (2) gij (gikgim — gugkm) {gjkgim — gjigkm) a -6 

= -48 (4) gij (2 9ij ) a" 6 = -48 (24) a" 6 

-4:4Rp Upa Rp Uti \Rp aK x = 44 (g ik gji - gugjk) (gi P gj q - gi q gj P ) (gk P gi q - gi q gi P ) a -6 
= 44 (4) (6) a' 6 

-80Rp Upcr Rp K p\R UKIJ \ = 80 (gikgji — gugjk) (9ik9 Pq — 9i q 9 P k) {9ji9 P q ~ 9j q 9 P i) a 6 
= 80 (Sgjigpg - g pq gji - gi } g vq + gi g g jp ) (gjig pg - gj q g p i) a~ 
= 80 (9 - 3 + 3 - 3) a~ 6 
= 80 (6) a~ 6 

Collecting the first set of terms we get 
4a-" 



9-7! 
3 



-35 • 6 3 + 42 • 72 - 42 • 72 + 208 • 24 - 192 • 24 + 48 • 24 - 44 • 24 - 80 • 6) 
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Now we continue with the second set of terms 

-12Tr (^^fi,,) = -12 Q Tr ( 7 „67 C( i7e/) R^R^R^ 

= Tr (1) 12 ^-^ (8) RnuabRvpbcRpfiac 

3 _ 6 
= - g (ff»fc5ji - 9il9jk) {gjlQpq ~ 9j P 9lq) (9 P k9i q ~ 9pi9kq) a Tr (1) 

3 _ 6 

= ~ 2 (^9ik9pq - 9ik9pq - 9ik9pq + 9iq9pk) {9pk9iq ~ 9pi9kq) 4a 

= - g (3-3 + 9-3) 4a" 6 
= -9 • 4a- 6 



6 

-6-R^po-Tr (f^Slpo-) = -—R^p^Ti (tabled) Rp V ah Rpa Cd 

= ~^RpvpaR[ 1 u Rpaab^Y (1) 

3 

= (fltaSji - (sVfljg - (#%>5z? ~ 9kq9l P ) 4a~ 

= -3(9-3) 4a~ 6 = -18 • 4cT 6 



-4i? M!/ Tr (ft^Ap) = — —R^Tr ( 7a67cci ) R^R^ 

= ^RpfRix P ab Rupab Tl (!) 



-R....R ah ' 
2 

1 _ 6 

2 (^) S'ij (,9ip9mq 9iq9mp) {9jp9mq 9jq9mp) 4a 



2 (9 - 3) 4a -6 = 12 • 4a 



-6 



5OTr (ny = ^OTr ( lablcd ) R^ ab R^ cd 
5 



8 -i?i?y CT Tr(l) 

= -^(6) (12) 4a- 6 
= -45 • 4a- 6 



30 



60Tr (E 3 ) = 60 (-^ R 3 Tt (1) 



60 (-^) 3 (6) 3 -4a- 6 
— (15-27) -4a- 6 



30Tr (M7 2 „) = 30 (-£) (-2) (-£) a^i^Tr 

= 1(5 -27) -4a- 6 
30 

301?£ 2 Tr (1) = — i? 3 Tr (1) 
u 16 w 

y) (6) 3 4a- 6 
= (15 • 27) • 4a~ 6 

5i? 2 £Tr(l) = --R 3 • 4 
4 

= -^(6) 3 4a" 6 
= - (10 • 27) • 4a~ 6 



(1) 



-2i? 2 ^Tr(l) = -2i^(- T )4 



-6 



R 
1 

-2(12) (~]4fl- 
36 • 4a" 6 



2i? 2 ^£Tr(l) = 2(12) {--) 4a- b 

= -36 • 4a' 6 

Collecting the second set of terms we get 

4a" 6 / 15 • 27 5 • 27 N 
9 + 18 - 12 + 45 + 15-27 + 10-27- 36 + 36 

OQU V Z Z j 



Thus the sum of all the terms in is zero. 
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